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We present Conedy, a performant scientific tool to numerically investigate dynamics on complex networks. 
Conedy allows to create networks and provides automatic code generation and compilation to ensure perfor- 
mant treatment of arbitrary node dynamics. Conedy can be interfaced via an internal script interpreter or 
via a Python module. 



Over the last decade, complex netw^orks research 
has contributed significantly to improve our un- 
derstanding of the dynamics of complex systems 
v^rhich are composed of mutually interacting com- 
ponents. The complicated relationship between 
network structure, dynamics of components, and 
emerging global behavior requires tailored tools 
that aid in unraveling new phenomena and in ver- 
ifying hypotheses and analytical results. Conedy 
(http://www.conedy.org) is a new scientific tool 
that allows to investigate dynamics on complex 
networks. In this paper, we give an overview of 
the functionality of Conedy and demonstrate its 
working procedure. 



INTRODUCTION 

A network is generally considered as a collection of 
nodes, which are connected by edges. The generality and 
simplicity of this notion has established the description of 
complex systems with networks as a standard approach 
in many scientific disciplines ^^"^ that are concerned with a 
large number of functionally similar and mutually inter- 
acting components. Examples range from sociology and 
quantitative finance via biology, earth and climate sci- 
ences to physics and the neurosciences. The availability 
of computers allows for numerical investigations of such 
complex networks on far larger scales than by a man- 
ual approach. Both, structural and functional aspects of 
natural and man-made networks have been investigated, 
often considering thousands or millions of nodes. 

One approach to describe complex networks is to pro- 
vide algorithms which generate model networks which 
are similar to complex networks in some aspect^^— . 
Apart from the merit of understanding structure and 
evolution of complex networks, such model networks 
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have inspired many studies of dynamics on networks, 
in which each node is associated with a dynamical 
system. These studies include investigations of dis- 
ease spreading2i2rii^ of neural network dynamics^"—, 
of noise-induced phenomena^^— and of synchronization 
phenomena32ri21. 

Edges in these networks typically indicate some kind 
of coupling, which is a weak disturbance of the node dy- 
namics. The question arises in which way the global net- 
work dynamics can be related to both the node dynam- 
ics and to the structure of the network. Can we observe 
global signals, which are similar to those of single ele- 
ments, or — even more interestingly — can we expect new 
emerging features which are not part of the dynamical 
spectrum of single elements? In order to untangle and 
understand the complicated relationship between net- 
work structure, global dynamics and node dynamics, a 
combined effort is necessary that includes theoretical ap- 
proaches and network^^Ti^ and time series analyses^"—. 
In this context, numerical studies play an important role 
in unraveling new phenomena, in verifying hypotheses 
and analytical results as well as in validating newly de- 
veloped network and time series analysis techniques. 

There are already software packages available that al- 
low to integrate different types of dynamical systems, but 
they either focus on neuron dynamics^^Ti^ or they do 
not consider networks^"—. Here, we present Conedy, a 
computational tool which is aimed at scientists investi- 
gating dynamics on complex networks. Conedy allows 
to build arbitrary networks and provides generators for 
many classic and popular networks, such as lattices, ran- 
dom, small world, and scale-free networks. Conedy can 
handle dynamical systems which are ordinary differen- 
tial equations, stochastic differential equations, iterated 
maps, and pulse-coupled oscillators. The dynamics of 
these systems can be assessed by a common interface, 
which easily allows to investigate the same network en- 
dowed with different dynamical systems as node dynam- 
ics. 

For ease of use, Conedy can be interfaced with 
an inbuilt script interpreter or via Python bindings 
I http : // www ■ python ■ org| ) . The inbuilt script inter- 
preter may be useful for distributed computing. It 



The following article has been accepted by Chaos. After it is published, it will be found at ,http : //link . aip . org/link/?cha 



2 



provides support for creating Condor dagman job files, 
thereby converting the iterations of loops into different 
jobs. Condor (|http : //www . cs .wise . edu/ condor|) is a 

job management system developed at the Computer Sci- 
ence Department of the University of Wisconsin. Python 
is a free, open source, cross-platform programming lan- 
guage and has become a viable alternative to Matlab and 
Octave recently^*. The standalone Python bindings of 
Conedy put a powerful scientific tool at the fingertips of 
the analyst, which can then be combined with an ever in- 
creasing number of other high-quality scientific packages 
for Python. 

Conedy is written in C-|— |- and licensed under the GNU 
license (|http : //www . gnu . org/ copylef t/gpl . html ) . 

Conedy is distributed as source code and as binary pack- 
ages for Windows and Linux (http: //www .cone dy . org, 
http://www.github.coiii/conedy/conedy|. Conedy is 
not intended to contain plotting routines or support 
for other forms of visualization, which are excellently 
handled by other free software packages. Before we give 
an overview of the functionality of Conedy we provide a 
short example to demonstrate its working procedure. 



A SHORT EXAMPLE 
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Figure 1. (Color online) Second variable of three coupled 
Rossler oscillators in the funnel regime. Top: coupling 
strength 0.055, bottom: coupling strength 0.2. 



In this example (written in Python), we simulate a 
small network of three coupled Rossler oscillators in the 
funnel regime, which display phase synchronization for a 
sufficiently large coupling strengt h^^'^" . The differential 
equations for an oscillator read: 



X = —ijjy — z 
y = Lox + ay 
z — b + z{x — c) 

The oscillators are non-identical (wi = 1.06, lu2 = 1.02, 
= 0.98, a = 0.22, b = 0.1, c = 8.5) and diffusively 
coupled via their y-components. The dynamics is inte- 
grated with the default integrator and the second dynam- 
ical variable of the oscillators {y) is sampled in periodic 
intervals. We add three Rossler oscillators to an empty 
network 



import conedy as co 






net 


= CO . network ( ) 






CO . 


set ( " roessler_a " , 0.22) 






CO . 


set ( " roessler_b " , 0.1) 






CO . 


set ( " roessler_c " , 8.5) 






rl 


= net . addNode ( CO . roessler () ) 






r2 


= net . addNode ( CO . roessler () ) 






r3 


= net . addNode ( CO . roessler () ) 






net 


. setParam (rl , " roessler_omega" , 


1 


06) 


net 


. setParam (r2 , " roessler_omega" , 


1 


02) 


net 


. setParam (r3 , " roessler_omega" , 





98) 


set 


random initial conditions. 






net 


. randomizeStates ( co . roessler () 








CO. uniform (-10.0, 10.0), 







CO. uniform (-5.0, 5.0), 
CO. uniform (-0.5, 1.5)) 



connect oscillator 1 bidirectionally to both 2 and 3 with 
coupling strength 0.2, 



net 


addEdge 


(rl , 


r2 , 


CO 


we ight edEdge 


(0 


2)) 


net 


addEdge 


(r2, 


rl , 


CO 


we ight edEdge 


(0 


2)) 


net 


addEdge 


(rl , 


rS , 


CO 


we ight edEdge 


(0 


2)) 


net 


addEdge 


(r3, 


rl , 


CO 


we ight edEdge 


(0 


2)) 



let 100 units of time pass to let transients die out, 



net . evolve (0 . , 100.0) 



select the time and each oscillator's second compo- 
nent (zero-based numbering) to be written to the file 
roessler. dat, 



net 


observeTime ("roessler 


dat " ) 




net 


observe (rl , "roessler 


dat " , 


CO . component (1)) 


net 


observe (r2 , "roessler 


dat " , 


CO . component (1)) 


net 


observe (rS , "roessler 


dat " , 


CO . component (1)) 



set the sampling interval to 0.01, and observe the system 
for 100 units of time. 



CO . set ( " samplingTime " , 0.01) 
net . evolve (100 . , 200.0) 



In Fig. [T]we show the second variable of the oscillators as 
calculated with the listed script and, for comparison, for 
a coupling strength 0.055. Phase synchronization can be 
observed for the larger coupling strength. In the follow- 
ing, we will describe the design of Conedy and exemplify 
its working procedure. 
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CHOOSING DYNAMICS 

As most networks considered in both theoretical and 
experimental studies are sparse, Conedy administers an 
internal adjacency list representation for the investigated 
network, which consists of nodes and edges of different 
type. Each node type is associated with a dynamics, 
which can be of different kind. Conedy can handle dy- 
namical systems of the following classes: 

ordinary differential equations: For numerical inte- 
gration of ordinary differential equations, Conedy 
uses algorithms of the GNU Scientific Library 
(GSL)^, which contains step-controlled algorithms 
like the widely used Runge-Kutta-Fehlberg method 
of orders 4 and 5. 

stochastic differential equations: Implementations 
of the Euler-Mayurama and the Milstein method 
are supplied^. 

pulse-coupled oscillators: Integration is handled by 
an event-based approaches.. A relaxed heap^ or 
a calendar queued can be chosen as priority queue 
at compile time. Interaction of oscillators is defined 
by phase response curves^. Conedy allows for de- 
layed interactions and Poissonian excitations. 

iterated maps: Implemented straight-forwardly. 

Conedy ships with a list of some predefined node dynam- 
ics such as Kuramoto-, Rossler-, and Lorenz-oscillators or 
different neuron niodela^^— . In case the desired dynam- 
ics is not in this list, Conedy supports automatic gener- 
ation of code using a small description file, which follows 
the syntax of Windows INI files. In this description file, 
the dynamics and an integrator has to be specified. In 
addition, control parameters for the dynamics can be de- 
fined. After compilation of the code a new node type is 
available in Conedy. At run time, different values can be 
assigned to the control parameters for nodes of this type. 

With the following example, we list and explain the 
description file, which would be required to add the dy- 
namics of an excitable Barkley unit^ to Conedy. The 
dynamics is often used in studies on pattern formation 
and is similar to that of a FitzHugh-Nagumo neuron. 
The following stochastic differential equation describes 
the dynamics which is driven by Gaussian white noise^^: 

dxo — yXoil — xo){xo ^——)dt + adW, 

k a 

dxi — {xq — Xi)dt. 
We start by defining a name for the node type: 

[gaussianBarkley] 

We choose the type of the dynamical system as stochastic 
differential equations. The desired stepping algorithm 
can then be specified in the Python script. 



type = sde 

Next, we specify the number of dynamical variables of 
the system, here the number of differential equations; 

dimension = 2 

We specify the number of control parameters for the dy- 
namics: 

parameter = 4 

Each control parameter is assigned a name and a default 
value. 

par amet ername 1 = a 
def ault value 1 = 0.75 
parametername2 = b 
def aultvalue2 = 0.01 
parameternameS = k 
def aultvalueS = 0.05 
par amet ername4 = sigma 
def aultvalue4 = 0.18 



Finally, we define the dynamics of the node type, i.e., 
the derivatives and as well as the noise terms 



dynami c s = 




dxdt [0] 


= l/k*x [0] * (1-x [0] ) * (x [0] - (x [1] +b)/a) ; 


dxdt [1] 


= x[0] - x[l] ; 


s[0] = 


s igma ; 


s[l] = 


0.0; 


dsdx [0] 


= 0.0; 


dsdx [1] 


= 0.0; 



As the noise is not multiplicative, the derivatives and 

can be omitted. The equations use C-syntax and 
may contain — in addition to standard C-constructs — 
mathematic functions as defined in math . h. 

While nodes are associated with dynamics, edges in- 
dicate some kind of coupling. The exact nature of this 
coupling depends on the node that the edge connects to. 
For example, for an iterated map the edge may signify 
the typical coupling for a coupled map lattice, while it 
may indicate diffusive coupling for ordinary differential 
equations. The coupling has to be defined alongside the 
internal dynamics of the node type. For this purpose the 
node is provided (by each edge) with two numbers, which 
are interpreted as state of the coupled node and as cou- 
pling weight. Edges of different type differ in the way 
they determine these two numbers. The weight may be 
stored per edge or only per edge type to save memory. 
The state may be — depending on the edge type — a vari- 
able of the coupled node (e.g. xq) or some transformation 
of a variable. 

To define the coupling, the macro f orEachEdge can be 
used, which loops over all adjacent edges, sets coupling 
weight and state to the values provided by this edge, and 
performs an instruction which is given as argument to 
the macro. To include diffusive coupling to the definition 
of the Barkley dynamics, we add the following line to the 
property dynamics of the description file: 
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f orEachEdge 

(dxdt[0] = dxdt [0] + weight *( state -x [0] ); ) 

The new node type is added to Conedy by depositing 
the description file into a directory which is monitored 
by Conedy. For sake of performance, Conedy has to be 
recompiled after each addition of a new node type. Under 
Linux this is handled automatically at the next import 
of Conedy to the Python interpreter. 

BUILDING NETWORKS 

In the following we explain the different approaches 
for network generation and manipulation available in 
Conedy. For the sake of simplicity, we assume that 
Conedy has been imported to Python in the standard 
way by 

import conedy as co 

and that a network has been declared by 

net = CO. network () 

At this point, the network net is empty and serves as 
a starting point for the various manipulation functions 
which allow to build the desired network step by step. 

Most of the functions in Conedy which manipulate 
networks accept node templates and edge templates as 
arguments. These are Python classes, which specify 
node/edge type and values of control parameters. To 
create a node template of a certain type using default 
control parameters, Conedy accepts the following syntax: 

nodeType () : Return a template of node type nodeType 
using default values for control parameters. 

Alternatively, it is possible to create a node template 
with user-defined values for control parameters: 

nodeType (pi, pM) : Return a template of type 

nodeType using pi, pM as values for the M 
control parameters the node type depends on. 

Templates for edges can be created in a similar way. 
Conedy supplies elementary functions, which add a single 
node or a single edge to the network: 

net . addNode (nt ) : Add a node to the network net ac- 
cording to the node template nt and return a 
unique identifying integer (the node number n) for 
the added node. 

net.addEdge(s, t, et) : Add an (directed) edge to 
the network net according to the edge template 
et connecting node s with node t. 

With the help of these elementary functions it is possi- 
ble to create arbitrary networks in Python. For building 
commonly used networks such as random networks or 
lattices, higher level functions are available, e.g.: 



net . randomMetwork(#, p, nt, et) : Add # nodes ac- 
cording to the node template nt to the network 
net. Connect each pair of newly added nodes with 
probability p by an edge according to the edge tem- 
plate et. Return the node number n of the first 
added node. Other added nodes have consecutive 
numbers starting with n + 1. 

net. line (iV, fc, nt , et) : Add N nodes according to 
the node template nt to the network net. Con- 
nect the newly added nodes in order to form an 
open chain in which nodes (except the ones near 
the boundaries) are connected to their fc nearest 
neighbors to each side. Return the node number n 
of the first added node. Other added nodes have 
consecutive numbers starting with n-l-1. 

In addition, small-world or random networks with a 
predefined degree distribution can be created using ma- 
nipulation functions that allow to add (or to replace 
existing edges by) random edges. Finally, it is possi- 
ble to create networks from adjacency lists or matrices 
stored in external files. To verify and investigate the 
created network, some standard network measures from 
network theory like clustering coefficient, average short- 
est path length and centralities are implemented. In ad- 
dition, Conedy provides export functions to transfer net- 
works to other network analysis tools such as igraph^ or 
NetworkXS. 



SETTING CONTROL PARAMETERS AND INITIAL 
CONDITIONS 

For networks which consist of nodes with one or a few 
sets of control parameters and initial conditions, it is con- 
venient to set these when defining the node templates. 
Additionally, it is possible to draw control parameters 
and initial conditions from random distributions (uni- 
form, Gaussian, exponential). Random number genera- 
tion for this (and for the stochastic integration) is han- 
dled by the GSL. Both algorithm and random seed can 
be controlled to allow for reproducible computations. 

net . randomizeParameter (param^/ame, ranDist) : 
Draws for every node — which is of the node type 
to which control parameter paramName belongs — a 
random value for paramName from the distribution 
ranDist. 

In a similar way, initial conditions of nodes in the network 
can be drawn from random distributions. 

net . randomizeStates (nt, ranDist{, ranDist }) : 
. Draw a random number for every dynamical vari- 
able of all nodes which match the node template 
nt. For every dynamical variable, a distribution 
ranDist has to be chosen. 
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CHOOSING OBSERVABLES AND STARTING 
NUMERICAL INTEGRATION 

Conedy writes the values of one or more selected ob- 
servables of the network (e.g., first dynamical variable 
of the second node, third dynamical variable of the first 
node, etc.) to one or several files. In addition, observ- 
ables such as mean network activity or Kuramoto's order 
parameter— can be selected. Conedy administers a list 
of observables for every file and orders the columns in 
these files in the same order as observe commands have 
been issued: 

net. observe (n, fileName, et): Add the state of 
node n — as returned by an edge according to tem- 
plate et — to the list of observables to be written 
to file fileName. 

The standard output is a matrix with rows corresponding 
to sampling times and columns corresponding to observ- 
ables of the network. Eventually, a call of the evolve 
function starts the integration, thereby sampling all reg- 
istered observables with a globally defined sampling in- 
terval: 

net. evolve (s, t) : Set the time to s and evolve the 
dynamics on the network net until the time t. 

Data can be written to whitespace-newline-separated 
text files, bzip-compressed text files, or to binary files. 



EXAMPLE #1: SMALL-WORLD NETWORK OF 
DIFFUSIVELY COUPLED STOCHASTIC BARKLEY 
UNITS 

As a first usage example, we consider the formation of 
wave patterns in a two-dimensional small-world network 
of diffusively coupled Barkley units driven by Gaussian 
white noise^i^. For this purpose we use the node type 
defined above. We begin with creating a node template 
and specify its initial conditions as the fixed point at 
{xo = 2/0 = 0.0). 



bark = co . gaus s i anBar kl ey ( ) 
bark . setState (0.0, 0.0) 



Next, we create a two-dimensional lattice of 512 x 512 
units, where each pair of units is connected if their Eu- 
clidean distance is smaller than or equal to 1.0 (nearest 
neighbor coupling) . The coupling strength for each edges 
is 3.84. The network is then rewired by replacing one per 
mill of edges with edges between randomly chosen source 
and target nodes. 



net . lattice (512 , 512, 
net . rewire (0 . 001) 



1.0, bark, 



. weightedEdge 



By applying a Milstein method to the resulting 524288- 
dimensional stochastic differential equation system, we 
let the time evolve to 20.0. This is sufficient to let tran- 
sients die out. 
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Figure 2. (Color online) Spiral waves in a small-world network 
of 512 X 512 diffusively coupled noisy Barkley units (values 
range from 0.0 (white) to 1.0 (black)). Snapshot of the spatial 
distribution of the first dynamical variable. 



net . evolve (0 . , 20.0) 



Next, we add the first dynamical variable of all nodes to 
the list of observables and specify the name of the file as 
waves . dat. 



net.observeAll ("waves.dat", co. component (0) ) 



Finally, we generate a snapshot of all registered observ- 
ables at the current time (20.0) of the network: 



net . snapshot () 



The data is then written to the file waves.dat and is 
shown in Fig. [S] At the time of the snapshot, spiral 
waves can be observed from the spatial distribution of 
the first dynamical variable of the units. 



EXAMPLE #2: RANDOM NETWORK OF 
PULSE-COUPLED OSCILLATORS 



As a second usage example we consider chaotic tran- 
sients to complete synchrony in a random network of 
pulse-coupled oscillators (PCOs)^'''^^. PCOs are de- 
scribed by their phase (j) G [0,1] with (p — I. When <j) 
reaches 1 the oscillator fires and (p is reset to 0. The in- 
teraction between PCOs is defined by the phase response 
(3.%¥)9tion A(0), i.e., the dependence of the phase change 
A on the phase (j> of an oscillator due to the firing of a 
connected oscillator. 

Whenever at some time t an oscillator receives a pulse, 
its phase is updated according to 



0(i+) = 0(i)-l-A(0). 
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Figure 3. Spontaneous transition to complete synchrony in a 
random network of 1000 pulse-coupled oscillators. Temporal 
evolution of Kuramoto's order parameter r{t). 

PERFORMANCE CONSIDERATIONS 

In scientific computing, it is often possible to find spe- 
cialized solutions, which perform well for certain condi- 
tions but either have a narrowed functionality or perform 
worse for other conditions. For numerical integration of 
complex network dynamics, different kinds of optimiza- 
tion can be achieved if the network is required to be in 
some way homogeneous. A more general approach will 
allow for different edge/node types and control parame- 
ters. 



Homogeneity of node and edge types 

A simulation tool for dynamics on complex networks 
that consist of different edge/node types (mixed network) 
can be realized by making use of virtual functions and 
jump tables to determine the edge types at run-time. 
This also offers the convenient possibility to change the 
type of edges without recompilation of the source code. 
At the downside, however, the use of virtual functions 
may be unnecessary slow for networks consisting of iden- 
tical types for all edges of a node (homogeneous net- 
works). In addition, the jump tables for virtual functions 
consume a considerable amount of memory. 

To enable an efficient handling of both situations (ho- 
mogeneous and mixed networks), Conedy follows two dif- 
ferent approaches to the internal representation of edges. 
First, Conedy allows for nodes that can connect to edges 
of arbitrary type (mixed networks) by making use of vir- 
tual functions. The second approach enables an efficient 
handling of homogeneous networks. For this purpose, 
Conedy allows for nodes, for which the connecting edge 
type is specified at compile time. Numerical integration 
can be achieved without virtual function calls using these 
nodes. To create a specialized node type, which connects 
to edges of a predefined type, we specify this edge type 
in the description file. For example, we complement the 
description file for the Barkley dynamics by: 



Here, we consider the case of a phase response function 
A(0) which is limited to the value that is needed for the 
excited oscillator to fire immediately^. 

A(0) = min{a + b<j),l- </>}, </> £ (0, 1). 

We assume that oscillators are not excitable at the time 
of firing (A(0) = A(l) = 0). The dynamics of these 
PCOs can be integrated into Conedy by the following 
description file. As above we declare name, integrator, 
dimension, and control parameters of a node: 

[pcoMirollo ] 
type = pco 
dimension = 1 
parameter = 2 
par amet ername 1 = a 
def aultvaluel = 0.01 
parametername2 = b 
def aultvalue2 = 0.02 

Instead of a differential equation, the property dynamics 
now contains the phase change delta in dependence on 
the phase (phase): 

dynami c s = 

if (phase == 0) { 
delta = 0; 

} 

else { 

delta = a + b*phase; 

if (delta + phase > 1.0) 

delta = 1.0 - phase; 

} 

Next, we create a random network of PCOs with N — 
1000 nodes of type pcoMirollo, where nodes are con- 
nected with a probability of 0.01 by unweighted edges. 

net . randomNetwork (1000 , 0.01, 
CO. pcoMirollo () , 
co.edgeO) 

Initial conditions are drawn from the uniform distribu- 
tion over the interval [0, 1]: 

net .randomizeStates( CO. pcoMirollo () , 
CO . uniform (0 . , 1.0)) 

We reset the default control parameters for the oscilla- 
tors: 

CO . set ( " pcoMirollo_a" , 0.015) 
CO . set ( " pcoMirollo_a" , 0.045) 

As observables we register the sampling time and Ku- 
ramoto's order parameter r{t) — 1/iV e^'^*'^"^*^ | for 
all N nodes. r{t) yields 1 for complete synchrony and 
for asynchronous firing. Both observables are to be 
written to the file order.dat, and we let the time evolve 
1000.0 time units. 

net . observeTime ( " order . dat " ) 

net . observePhaseCoherence ("order. dat") 

net . evolve (0 . , 1000.0) 

Fig. 13] shows the temporal evolution of the order pa- 
rameter r{t), which exhibits a spontaneous transition to 
synchrony. 



staticEdgeType=component_weightedEdge 



With this example, nodes of type gaussianBarkley are 
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Figure 4. (Color online) Execution times in dependence on 
the mean degree of random networks of 5000 Rossler oscil- 
lators with eigenfrequency 0.89. Values were obtained using 
nodes allowing for inhomogeneous node/edge types (blue dot- 
ted line) and using nodes which force connecting edge types 
at compile time (black solid line) on a PC with 1.5 GFLOPS 
(2.4 GHz). Integration method: Runge-Kutta-Fehlberg with 
step size control and an absolute error of 10~^. The dynamics 
of the networks was integrated for 200.0 time units. (Lines 
are for eye-guidance only). 



connected by edges of type component_weightedEdge 
only. 

With the second approach the memory consumption 
of unweighted edges in homogeneous networks can be re- 
duced by a factor of four^^. To compare both approaches, 
we numerically integrate the dynamics of a random net- 
work of 5000 Rossler oscillators. In Fig. |4] we show the 
respective execution times in dependence on the mean 
degree: the higher the mean degree, the greater the per- 
formance benefit of the approach specialized for homoge- 
neous networks. 



Homogeneous node parameters 

For networks consisting of a large number of identical 
nodes with identical control parameters, it is not desir- 
able to reserve memory for the parameters of each node. 
However, a specialized implementation for fixed parame- 
ters will narrow functionality. To provide a reasonable 
implementation for both homogeneous and inhomoge- 
neous node parameters, Conedy uses parameter sheets, 
which are stored in a lookup table. Each node reserves a 
32 bit integer for the number of the sheet which contains 
the node's parameters. Nodes with identical parameters 
share the same sheet number which reduces memory con- 
sumption. When parameters of nodes are drawn from a 
random distribution, each node is assigned its own pa- 
rameter sheet. Note, that in this situation, the chosen 
implementation will lead to a slightly increased memory 
consumption as compared to a direct storage of control 
parameters as part of the nodes. 



Event-driven integration of pulse-coupled oscillators 

Conedy handles the numerical integration of pulse- 
coupled oscillators (PCOs) via an event-driven 
approach^, which makes use of a priority queue. 
This queue contains the PCOs ordered by the time of 
their next firing. When an oscillators fires, it is removed 
from the top of the queue and reinserted at the bottom. 
In addition, the priority (firing time) of all connected 
oscillators is decreased/increased. Depending on the 
implementation of the queue, these four operations 
(remove, insert, increase priority, decrease priority) may 
scale differently with the size of the queue. In principle, 
with the use of a calendar queued, constant scaling for 
all four operations is possible. However, this queue has 
a comparably large overhead and only performs well, 
if events are distributed homogeneously in time, which 
corresponds to asynchronous oscillators. Alternatively, 
heap structures can be used which are equally well suited 
for arbitrary event distributions. They lack, however, an 
efficient way to decrease the priority of events (which is 
necessary to handle negative phase responses). 

In Conedy, both a relaxed heap and a calendar queue 
can be chosen at compile time. Compared to a clock- 
driven approach (e.g., a Runge-Kutta scheme), the inte- 
gration is exact (as far as admitted by double precision) 
and does not require the choice of a step size. To compare 
an event-driven with a clock-driven integrator we mea- 
sure execution times for random networks of either 10000 
excitatory or 10000 inhibitory pulse-coupled integrate- 
and-fire neurons- In Fig. [S]we show execution times 
in dependence on the mean degree-^. We chose control 
parameters, such that neurons fire asynchronously. The 
relaxed heap outperforms the calendar queue for exci- 
tatory oscillators, while the calendar queue seems bet- 
ter suited for inhibitory oscillators. For large mean de- 
grees, the event-driven integration is in both cases out- 
performed by a clock-driven approach, in which we check 
for firing oscillators after every integration step. This is 
because in the phase description, excitations of integrate- 
and-fire neurons are computationally more costly as they 
require the evaluation of exponentials. 



CONCLUSION 

Conedy is a convenient, easy-to-use scientific tool for 
people interested in investigating dynamics on complex 
networks. It allows to create networks via elementary 
functions (which add single nodes and edges), via higher 
level functions for often investigated networks (e.g., lat- 
tices, random, small world, and scale- free networks), or 
from external data (in the form of adjacency matrices or 
lists). The size of a network is limited by the available 
computer memory only. Conedy provides automatic code 
generation and compilation to ensure performant treat- 
ment of arbitrary node dynamics. We here have concen- 
trated on the Python interface of Conedy and on usage 
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Figure 5. (Color online) Execution times per firing neuron 
T in dependence on the mean degree for random networks of 
10000 integrate-and-fire neurons. Values were obtained us- 
ing event-driven approaches implemented in Conedy with a 
calendar queue ("cq") blue), with a relaxed heap ("rh", red) 
and with a clock-driven approach ( "cd" , black; Runge-Kutta- 
Fehlberg with step size control and an absolute error of 10~^). 
The calculation was done for both inhibitory and excitatory 
neurons. While the relaxed heap implementation is much 
faster for excitatory neurons ("rh, excit."), both types of neu- 
rons are integrated equally fast with the calendar queue and 
the clock-driven integrator. Tests were performed on a PC 
with 1.5 GFLOPS (2.4 GHz). (Lines are for eye-guidance 
only). 



examples. 

As algorithms for numerical integration have become 
complex and at the same time very reliable, it is prudent 
to hide much of their complexity behind frameworks to 
clear the sight for the system under study. Conedy com- 
bines several integration schemes that are summarized by 
a common interface, which easily allows to investigate the 
same network endowed with different kinds of dynamical 
systems as node dynamics. As with similar scientific soft- 
ware packagesAS-i^^ , it is advisable that the potential user 
has a thorough background in the numerical treatment 
of dynamical systems which cannot be provided by this 
paper, but rather by textbooka^i^S. The user will still 
need considerable experience in order to avoid spurious 
interpretation of the results, especially for large networks. 

The last years have seen an extraordinary success of 
complex network and their applications in diverse disci- 
plines. We hope that with Conedy we can contribute to 
this rapidly evolving field. 
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